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This  document  contains  the  slide  presentation  entitled  "High-Frequency 
Acoustic  Backscatter  From  the  Sea  Surface,*  given  at  the  107th  meeting  of  the 
Acoustic  Society  of  America  on  8  May  1984  in  Norfolk,  Virginia.^ 

"A  high  resolution  scattering  experiment  was  conducted  in  the  shallow  waters 
of  the  North  Atlantic.  A  narrow  beam  parametric  array,  which  was  rotatable  in 
both  azimuth  and  elevation,  was  utilized  as  a  broadband  high-frequency  acoustic 


DO  l  , *N*7J  1473  EDITION  of  I  NOV  «t  IS  OBSOLETE 
S/d  0  102-014-  it 0  1 


20.  Continued: 


projector.  Acoustic  surface  scattering  data  were  obtained  at  normal  incidence 
and  low  grazing  angles  (less  than  10°)  as  a  function  of  acoustic  transmit 
frequency  and  sea  state  conditions.  Meteorologic  and  oceanographic  data  were 
obtained  in  concert  with  the  acoustic  measurements  and  included  wind  speed  and 
direction,  ocean  surface  wave  spectra  and  currents,  and  ocean  sound  speed. 
Surface  backscatter ing  strength,  Doppler  spectra  (shift  and  spread),  and 
envelope  statistics  were  some  of  the  measured  parameters.  It  will  be  shown  that 
the  Doppler  spectra  are  approximately  Gaussian  and  the  spectral  shift  could 
be  predicted  from  Bragg  diffraction  theory  modified  by  the  induced  Doppler 
due  to  surface  currents.  The  measured  Doppler  spreads,  however,  were  greater 
than  those  predicted  by  composite  roughness  theory^  At  low  grazing  angles,  the 
surface  backscatter ing  strength  showed  a  lack  of  s£hqng  frequency  dependence. 
Yet,  at  normal  incidence  under  nearly  constant  sea  surface  conditions,  the 
surface  loss  varied  20  dB  as- the  frequency  changed  from'S.  to  80  kHz,  which 
disagrees  with  physical  optics  in  the  high  frequency  limitV 
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HIGH-FREQUENCY  ACOUSTIC  BACKSCATTER 
FROM  THE  SEA  SURFACE 


INTRODUCTION 

A  high  frequency  acoustic  boundary  reverberation  experiment  has  been 
conducted  In  the  shallow  waters  of  the  North  Atlantic.  The  experiment  was 
jointly  performed  among  several  Naval  centers  and  academic  Institutions,  In¬ 
cluding  NORDA;  FWG,  Kiel,  Germany;  the  University  of  Rhode  Island;  City 
College  of  New  York;  and  NUSC. 

In  a  shallow  water  environment,  propagation  under  refractive  or  nonre- 
fractlve  conditions  generally  Involves  reflection  and  scattering  from  the 
ocean  boundaries.  Interaction  with  the  ocean  boundaries  generally  occurs  at 
low  grazing  angles;  and,  because  of  the  scarce  and  often  conflicting  back- 
scatter  data  that  have  been  obtained  previously,  a  primary  objective  of  the 
experiment  was  to  determine  the  statistical  characteristics  of  the  surface 
reverberation.  An  attempt  was  also  made  to  measure  the  normal  incidence 
backscatter  and  Interpret  these  data  In  terms  of  the  ocean  wave  parameters 
and.  If  present,  near-surface  bubble  layer  thickness  and  volume  scattering 
strength.  Bottom  backscatterlng  measurements  made  during  this  experiment  are 
described  In  a  companion  paper  [1]. 

Presented  In  this  paper  are  some  of  the  envelope  statistics,  scattering 
strengths  as  a  function  of  frequency,  spectral  characteristics  of  the  sur¬ 
face  reverberation,  and  surface  wave  parameters.  As  will  be  seen.  It  was 
equally  Important  to  measure  the  oceanographic  conditions  and  statistics  of 
the  ocean  surface  for  subsequent  modeling  of  the  acoustic  data. 


Figure  1 .  Sea  Surface  Backscatter  Experiment 
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ACOUSTIC  ANO  OCEANOGRAPHIC  MEASUREMENTS 

A  parametric  array  (figure  1)  was  mounted  on  the  top  of  a  7  m  high  plat¬ 
form  and  was  projected  at  the  sea  surface  at  angles  of  normal  incidence  and 
grazing  angles  down  to  2*.  The  projector  was  also  rotated  in  azimuth  to 
obtain  acoustic  backscatter  for  different  orientations  to  the  surface  wave 
directionality.  A  wave  rider  buoy  was  used  to  collect  ocean  wave  statistics. 
Including  directionality  throughout  the  experiment.  Wave  statistics  were 
also  obtained  at  normal  Incidence  backscatter  by  an  acoustic  Inversion  tech¬ 
nique  developed  previously  [2].  Both  wave  measuring  techniques  were  limited 
to  resolving  surface  wavelengths  greater  than  approximately  1  m.  The  near 
surface  current  array  was  used  to  obtain  estimates  of  the  surface  current 
magnitude  and  direction.  For  In-sltu  ray  arrival  determination  and  subse¬ 
quent  propagation  loss  prediction,  conductivity,  temperature,  and  depth 
measurements  were  performed  periodically  during  the  acoustic  experimentation. 
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PULSE:  1  A3  ENSEMBLE:  1  to  400  PULSE  LENGTH:  0.4  imoe 
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Figure  2.  Normal  Incidence  Surface  Backscatter  Envelope 
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NORMAL  INCIDENCE  SURFACE  BACKSCATTER 

Acoustic  backscatter  from  the  sea  surface  was  obtained  at  normal  inci¬ 
dence  using  a  parametric  array,  which  also  functioned  as  a  conventional 
receiver.  The  insonified  area  on  the  surface  was  on  the  order  of  one  square 
meter.  The  transmitted  pulses  were  of  Identical  shape  and  phase  having  the 
same  zero  time  reference  for  each.  The  upper  three  plots  in  figure  2  show 
the  envelope  time  history  of  three  Individual  received  pulses  for  a  frequen¬ 
cy  of  20  kHz  and  a  0.4  ms  pulse  length.  Fluctuations  In  received  amplitude 
on  a  pulse-to-pulse  basis  are  a  result  of  the  time  variation  In  the  average 
number  of  specular  points  on  the  surface  times  the  average  curvature  of 
these  points  as  predicted  from  either  geometrical  or  physical  optics  theory. 
The  bottom  figure  Is  the  ensemble  average  of  400  pulses.  The  time  dispersion 
In  the  ensemble  is  approximately  3  ms,  which  Is  directly  related  to  the 
vertical  height  variation  in  the  surface  wave  motion. 
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In  order  to  use  acoustic  Inversion  techniques  to  measure  the  ocean  wave 
parameters,  near-surface  bubble  scattering  layer  thickness,  and  volume 
scattering  strength,  0.4  ms  pulses  were  transmitted  and  ensemble  averaged. 
The  0.4  ms  pulses  would  resolve  bubble  layer  thickness  to  0.75  m.  A  thresh¬ 
old  level  of  -20  dB  was  set  to  determine  the  onset  of  surface  backscatter 
and  400  pulses  were  ensemble  averaged.  The  result  Is  shown  in  figure  3  for  a 
transmitted  frequency  of  20  kHz.  There  Is  no  clear  evidence  of  a  persistent 
near-surface  bubble  scattering  layer,  which  would  show  as  a  precursor  to  the 
onset  of  the  surface  backscatter.  The  minimum  detectable  volume  scattering 
strength  was  -47  dB. 
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Figure  4.  Surface  Scattering  Coefficient  Vs 
Rayleigh  Parameter  at  Normal  Incidence 
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The  normal  incidence  surface  scattering  coefficient  was  calculated  using 
a  physical  optics  approach  [3]  and  Is  shown  In  figure  4  as  a  function  of  the 
Rayleigh  parameters  for  several  acoustic  frequencies  and  two  sea  state  con¬ 
ditions.  At  these  large  Rayleigh  parameters,  the  acoustic  intensity  Is  Inde¬ 
pendent  of  both  frequency  and  root-mean-square  wave  height  and  is  solely  a 
function  of  the  mean  squared  slope  and  Insonlfled  area.  It  can  be  seen  that 
there  Is  considerable  discrepancy  between  the  measured  values  and  the  theo¬ 
retical  prediction.  The  coefficient  of  variation  for  some  of  these  data  sets 
was  greater  than  52. 5%,  which  Indicates  that  the  Insonlfled  area  was  smaller 
than  the  surface  correlation  length. 
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ELEVATION  ANGLE:  »0  VARIANCE  946.7  em2  RMS  HEIGHT  30.8  om 
PULSES  1  to  400  MEAN  0.0  em  WIND  SPEED  19.4  kt 


0  16  32  48  64  80  96  112  128  144  160 


TIME  (aec) 

Figure  5.  Wave  Height  Vs.  Time 
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ACOUSTIC  WAVE  HEIGHT  MEASUREMENTS 


The  time  varying  wave  height  (figure  5)  was  calculated  from  the  travel 
time  associated  with  the  threshold  detected  surface  backscatter  and  the  mean 
sound  speed.  The  surface  height  was  sampled  at  a  rate  of  2.5  times  a  second, 
which  was  the  repetition  rate  of  the  transmitted  pulses.  The  standard  devia¬ 
tion  of  the  time  varying  wave  height  was  30.8  cm. 


The  probability  density  function  (PDF)  for  the  previous  wave  height  time 
history  Is  shown  In  figure  6  as  a  histogram  normalized  by  Its  area.  The 
solid  curve  Is  Gaussian  with  zero  mean  and  area  one.  Note  that  the  measured 
PDF  Is  slightly  skewed  toward  negative  heights  reflecting  the  fact  that 
waves  spend  more  time  as  troughs  than  peaks. 


MAGNITUDE  (cm2  mc) 


The  surface  wave  spectrum  was  calculated  from  the  wave  height  time 
series  derived  from  the  thresholded  acoustic  signals.  Comparisons  of  the 
mean  squared  wave  height  obtained  from  this  acoustic  technique  and  the  wave 
rider  buoy  showed  very  good  agreement.  This  verification  was  important  be¬ 
cause  the  wave  rider  buoy  data  were  obtained  on  an  hourly  basis  throughout 
the  experiment.  The  peak  of  the  spectral  energy  shown  in  figure  7  occurs  at 
a  frequency  of  0.2  Hz  and  the  root  mean  square  height  was  30.8  cm.  A  further 
comparison  with  a  Pierson  Moskowltz  [4]  spectrum  shows  that  for  this  wind 
speed,  the  sea  was  not  fully  developed  since  the  peak  frequency  should  be  at 
0.14  Hz  and  the  corresponding  wave  height  should  be  53  cm.  This  points  out 
the  necessity  of  measuring  surface  statistics  during  the  acoustic  measure¬ 
ments  and  not  relying  on  local  wind  speeds  as  an  Indirect  Indicator  of  sea 
conditions. 


ENVELOPE  (dB//1  volt) 
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SURFACE  REVERBERATION  AT  A  LOW  GRAZING  ANGLE 

A  low  grazing  angle  example  of  surface  reverberation  at  20  kHz  Is  repre¬ 
sentative  of  the  time  dispersion  and  amplitude  variability  commonly  observed 
for  backscatter  from  the  sea  surface.  The  top  three  plots  In  figure  8  are 
received  levels  from  consecutive  100  ms  pulses  obtained  with  a  repetition 
period  of  3.5  s.  The  bottom  plot  is  an  ensemble  average  of  the  first  four 
pulses  In  the  sequence.  As  expected,  the  averaging  process  reduces  signifi¬ 
cantly  the  amplitudlnal  Interference  effects  seen  in  the  individual  returns. 
The  9°  surface  grazing  angle,  which  corresponds  to  the  MRA  of  the  parametric 
source,  occurs  at  a  time  of  approximately  230  ms  In  this  data  set.  Because 
of  the  existing  sound  speed  conditions,  the  MRA  was  downwardly  refracted  by 
4#  before  Intersecting  the  sea  surface. 
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GRAZING  ANQLK:  ••  PULSKB:  1  to  250  TIME:  228  mooe 


Figure  9.  PDF  of  8ackscattered  Levels  for  a 


Fixed  Travel  Time 
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The  PDF  of  the  variations  In  the  envelope  of  the  backscattered  levels 
have  been  computed  from  the  surface  reverberation  measurements.  An  example 
of  a  9°  grazing  angle  (figure  9)  was  obtained  at  20  kHz  from  250  levels 
associated  with  the  travel  time  of  the  MRA.  Theoretically,  a  Rayleigh  dens¬ 
ity  is  predicted  for  the  envelope  fluctuations  of  incoherently  scattered 
energy.  The  measured  coefficient  of  variation  for  this  example  is  58%,  which 
compares  favorably  with  the  52%  expected  for  the  Rayleigh  density  function. 
Additional  statistical  tests  are  being  conducted  on  these  data.  PDFs  for 
single  ping  statistics  have  also  been  computed. 
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Surface  reverberation  was  measured  as  a  function  of  frequency  in  5  kHz 
intervals  from  5  to  20  kHz.  The  grazing  angle  was  held  constant  at  9°.  Dur¬ 
ing  the  measurements,  the  significant  wave  height  varied  from  0.7  to  1.0  m. 
It  can  be  seen  from  figure  10  that  the  scattering  strength  has  less  than  2.5 
dB  variation  when  the  frequency  Is  varied  over  two  octaves.  The  average 
scattering  strength  is  approximately  -30  dB.  A  comparison  with  Chapman  and 
Harris'  semiempirlcal  equation  for  scattering  strength  would  show  7  dB  vari¬ 
ation  over  the  same  frequency  range. 
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DOPPLER  CHARACTERISTICS  OF 
SEA  SURFACE  REVERBERATION  SPECTRA 

Shown  in  figure  11  are  sea  surface  reverberation  spectra  for  a  20  kHz 
transmit  frequency,  a  100  ms  pulse  length,  and  a  9°  grazing  angle.  The  zero 
on  the  frequency  axis  corresponds  to  the  energy  at  the  transmit  frequency. 
The  upper  three  spectra  demonstrate  the  variability  in  the  spectral  content 
for  three  consecutive  pulses,  each  implying  a  slightly  different  velocity 
distribution  of  sea  surface  scatters.  The  bottom  spectrum  is  an  ensemble 
average  of  85  Individual  spectra,  and  exhibits  an  approximately  Gaussian 
shape.  Note  that  the  average  spectral  peak  of  the  surface  reverberation  is 
downshifted  12  Hz  from  the  transmitted  frequency  and  the  energy  is  spread. in 
frequency  over  a  bandwidth  of  approximately  24  Hz. 
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The  Doppler  shift  of  the  average  spectral  energy  can  be  predicted  from 
resonance  theory  backscatter  as  suggested  by  Bass  and  Fuks  [5,6]  and 
Kur'yanov  [7].  In  figure  12,  we  have  taken  Into  account  the  directions  of 
the  surface  wave  propagation  and  surface  current  relative  to  the  orientation 
of  the  source/receiver.  The  Doppler  shift  Is  predicted  from  the  surface 
wavelength  component  that  resonates  with  the  acoustic  wavelength  for  a 
given  grazing  angle  and  from  the  surface  current  component.  It's  Interesting 
that  the  Ooppler  shift  Is  Independent  of  surface  amplitude,  however,  this 
must  be  taken  Into  account  when  predicting  the  Doppler  spread. 


DOPPLER  SHIFT  (Hz) 
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Doppler  shift,  measured  from  the  surface  reverberation  spectra.  Is 
plotted  versus  transmit  frequency  In  figure  13.  The  theoretical  shift  (cal¬ 
culated  from  the  previous  equation,  and  the  measured  geometric,  oceanogra¬ 
phic,  and  acoustic  parameters)  Is  plotted  as  the  solid  line.  It  Is  evident 
that  the  measured  shifts  are  well  described  by  the  theoretical  equation. 
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DOPPLER  SPREAD  (Hz) 
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Figure  14.  Surface  Reverberation  Ooppler  Spread:  Theory  and  Measurement 
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Doppler  spread  (shown  In  figure  14)  was  measured  from  the  Doppler  rever¬ 
beration  spectra  between  the  10  log  1/E  points  below  the  spectral  maximum. 
The  theoretical  solid  curve  was  derived  from  an  expression  for  spectral 
spread  by  Bass  and  Fuks  [6].  The  predicted  spectral  spread  of  the  sound 
backscattered  from  the  sea  surface  is  Inversely  proportional  to  acoustic 
wavelength  and  orbital  period  for  the  maximum  surface  wave  energy  density. 
As  the  signal  frequency  Increases,  the  observed  spectral  spread  Increases  as 
predicted  by  the  expression.  The  theoretical  spectral  width,  however,  is 
only  approximately  40-60X  of  our  measured  width.  Bass  and  Fuks  compared 
their  spread  equation  with  measured  acoustic  data  and  observed  the  same 
result. 
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SUMMARY 


•  DOPPLER  SHIFTS  PREDICTED  BY  BRAGG  DIFFRACTION  THEORY 
WITH  SURFACE  CURRENT 


•  DOPPLER  SPECTRA  WERE  GAUSSIAN 


•  DOPPLER  SPREADS  WERE  GREATER  THAN  PREDICTED  BY 
COMPOSITE  ROUGHNESS  THEORY 


•  NORMAL  INCIDENCE  BACKSCATTER  DATA  DISAGREES  WITH 
PHYSICAL  OPTICS  HIGH  FREQUENCY  LIMIT  THEORY 


•  LACK  OF  STRONG  FREQUENCY  DEPENDENCE  (5-20  kHz)  OF 
SCATTERING  STRENGTH  AT  LOW  GRAZING  ANGLES 


•  ACOUSTIC  INVERSION  TECHNIQUES  SUCCESSFULLY  OBTAINED 
OCEAN  WAVE  PARAMETERS  AND  STATISTICS 


Figure  15.  Surface  Reverberation  Summary  Highlights 
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SUMMARY  AND  CONCLUSIONS 

The  experiment  has  verified  that  the  Doppler  shift  of  ensemble  averaged 
surface  reverberation  spectra  can  be  predicted  from  resonance  backscatter 
theory  and  that  the  ensemble  average  exhibits  an  approximately  Gaussian 
shape.  The  major  conflicts  between  theory  and  measured  data  occur  in  the 
prediction  of  (1)  Doppler  spread  at  low  grazing  angles  and  (2)  the  absolute 
levels  of  backscattered  intensity  as  exemplified  by  the  normal  incidence 
data.  The  high  values  of  scattering  strength  and  lack  of  strong  frequency 
dependence  cannot  be  predicted  from  current  backscatter  composite  models. 
This  can  also  be  seen  from  the  Doppler  spread  predictions  that  underestimate 
the  bandwidth  of  the  reverberation  Doppler  spectra.  As  an  aside,  it  is  en¬ 
couraging  to  see  that  acoustic  inversion  techniques  can  be  used  to  obtain 
important  statistical  characteristics  of  ocean  waves. 
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